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1. Introduction 



The high degree of isotropy of the cosmic microwave background radiation is usually 
taken as evidence that the late (matter-dominated) universe is isotropic about our 
worldline. When combined with the Copernican principle that we do not occupy a 
privileged position, this leads to isotropy about all worldlines, and thus to a Friedmann- 
Robertson- Walker (frw) universe. (The proof of this background-radiation argument 
is the Ehlers-Geren-Sachs theorem M, whose perturbed version is proved in 0.) 

If we suspend the Copernican assumption in favour of a directly observational ap- 
proach, as set out in detail in ||, then it turns out that isotropy of the background 
radiation is insufficient to force isotropy onto the spacetime geometry. However, it has 
been shown || that isotropy of galactic number counts and area distances (equivalently 
luminosity distances), together with the vanishing of cosmological proper motions and 
image distortion, does force isotropy of the universe about our worldline. The uni- 
verse is then a Lemaitre-Tolman-Bondi (ltb) model ||, including the special case of 
homogeneous frw models. 

Current galactic observations are nowhere near the accuracy and completeness of 
the background radiation observations, but they are not inconsistent with isotropy. 
However, observations of the galactic distribution do not imply homogeneity 0.0; 
indeed, they sometimes appear to indicate that inhomogeneities may not be smoothed 
out on larger and larger averaging scales. Therefore, motivated by observations, and 
leaving aside the Copernican principle, we assume isotropic galactic observations which 
are not necessarily (i.e. without further observational evidence) homogeneous. Obser- 
vational characterisations of homogeneity are then needed in order to test the basis on 
which the standard frw cosmologies rests. In it was shown that an ltb universe 
is homogeneous if and only if the observed galactic area distances and number counts 
take precisely the frw form (as functions of redshift). In Section 4, we derive a further 
characterisation of homogeneity within the observational framework, and show that 
for the parabolic case, the characterisation of || may be weakened: the single (area 
distance, number count) relation must take the Einstein-de Sitter form. 

The 3 + 1 analysis of ltb spacetimes ||, based on spatial hypersurfaces {t = 
const.}, is effective for finding and analysing exact solutions, but is not directly suited 
to observational cosmology. Since observations lead to data not on {t = const.}, but on 
the past light cone of the observer, a 2 + 2 approach is directly adapted to the analysis 
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of ltb universes from the viewpoint of observational cosmology. Furthermore, it is 
not possible simply to transform the 3 + 1 ltb solutions to observational coordinates, 
since such a transformation requires solution of the null geodesic equation, which is 
in general only possible numerically. In observational coordinates, the null geodesies 
are known directly, but the field equations are considerably more complicated, and the 
analogue of the 3 + 1 exact solutions is not known. In Section 2, we briefly review the 
observational coordinate approach. 

This paper complements and corrects ||, where it is indicated how in principle the 
field equations can be integrated using the redshifts, number counts and area distances 
as data. An important technical mistake was made in || (see Section 4), which severely 
over-restricts the solutions obtained by that procedure. Here, in Section 3, we give an 
explicit reduction of the field equations in 2 + 2 observational coordinates to a single 
first-order ordinary differential equation (ode) on the past light cone. The variables 
in the ODE are the area distance and a number count function. This corrects the 
integration procedure of ||, and is parallel to the 3 + 1 case of an explicit general 
solution. 

An application of this reduction is the construction in Section 5 of a class of explicit 
exact solutions representing an observationally-based inhomogeneous generalisation of 
flat frw universes. The solutions approach the Einstein-de Sitter solution near the 
observer, so that they are locally homogeneous (for small redshift). At greater redshifts, 
the solutions deviate appreciably from the Einstein-de Sitter solution: the matter 
distribution is locally homogeneous, but develops inhomogeneity at larger redshifts. 

2. Metric and observations in observational coordi- 
nates 

In 3 + 1 coordinates based on spatial hypersurfaces {t = const.}, the metric is || 



ds 2 = -dt 2 + 



dR(t,r) 
dr 



2 



dr 



2 



1 - kf(r) 2 



+ R(t,r) 2 dtt 2 (1) 



where k = 0,1, —1 corresponds, respectively, to parabolic, elliptic and hyperbolic in- 
trinsic geometry of the {t = const.} hypersurfaces, and dQ 2 is the metric of the unit 
2-sphere. The coordinates are co moving with the dust four- velocity = 5^ t - If r is 
chosen as a proper area coordinate, then the frw case is given by 

R(t,r) = ra(t) , f(r) = r 
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In the general case, the solution R(t, r) of Einstein's equations may be given explicitly 
||. However these coordinates and solution are not adapted to observations, which are 
made not on a {t = const.} surface, but rather on our past light cone. (See fH . |§ . || 
for further discussion of the crucial implications of this point.) Light emitted from 
a source at radial distance r E at cosmic time t# propagates along the null geodesic 
\/l — kpdt = —(dR/dr)dr before being observed at r = at time to > £e- The 
explicit analytic form t = t(r) of the null geodesic cannot in general be determined. 
If we use observational coordinates based on the past light cones of the observer's 
worldline, then the observations are given in the most simple and direct form. The 
price paid for this is that the field equations can no longer be solved explicitly. However, 
as we shall show, they may be reduced to a single first-order ode. 

In observational coordinates x M = (w, y, 9, 0), the metric takes the 2 + 2 form || 

ds 2 = -A(w, yfdw 2 + 2A(w, y)B(w, y)dwdy + C(w, yfdfl 2 (2) 

where {w = const.} are the past light cones along {y = 0}, and y is a comoving radial 
distance parameter down the light rays (null geodesies) {w = const., (9, 0) = const.}. 
Regularity on {y = 0} imposes on the metric the central limiting conditions ||: as 
y^O 

A(w,y) = A(w,0) + O(y) A(w,0)^0 

B(w,y) = B(w,0) + O(y) B(w,0)^0 (3) 
C(w,y) = B(w,0)y + O(y 2 ) 

The nontrivial coordinate freedom is a scaling of w and y: 

~ ~ /dw dy 

w ->• w = w{w) , y^y = y{y) (— t t -j-) ( 4 ) 

The first corresponds to a freedom to choose w as any time parameter (e.g. proper 
time) along {y = 0}. The second corresponds to a freedom to choose y as any null 
distance parameter (e.g. affine parameter, redshift) on an initial light cone {w = Wq}, 
after which y is dragged along by the matter flow. These scalings imply a scaling of 
the metric functions A, B: 

t dw . „ ~ dy „ 
A A = —A , B -> B = -^B 
dw dy 

The metric function C is the observer area distance (or angular diameter distance) 
on each past light cone, since dSs = C 2 dQE, where dSs is the cross sectional area of 
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the source and (EIe the solid angle it subtends at the central observer {y = 0} 0. If 
(1a denotes the area distance that is measured here-and-now, then 



dA(y) =C(w ,y) (5) 

where y may be chosen as the redshift, so that observations determine C on {w = w^}. 
By contrast, in the 3 + 1 approach, d,A_{r) = R(t(r),r), where R is known from the 
explicit exact solution, but r is not known and neither is the explicit form t(r) of the 
past light cone of here-and-now. The luminosity distance is [[| 

d L = (l + zfd A (6) 
where the redshift z of the source relative to the central observer is given by M: 

A{w ,y) 

The metric function B determines the deviation of y from the affine parameter v 

Two four-vectors are defined intrinsically by the matter and light rays: the dust velocity 
and photon wave vectors 

u» = A~ l 5^ w , P = (AB)- l 5» y (9) 

The number of galactic sources counted by the central observer out to a distance y 
is [see 0, which corrects the error in 0, equation (19)] 

N(y) = 4tt f n(w , y)B(w , y)C(w , yfdy (10) 
Jo 

on using © and (|). In (0), select ion effects are ignored, the freedom (||) has been 
used to set A(w , 0) = 1, and n is the number density of sources, so that the matter 
density is 

p = 47T77TO 

where Anm is the average galactic mass. 
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The rate of expansion of the dust is 3H = so that, by (Q), 

where a dot denotes d/dw. For the central observer, relative to whom the recession of 
other worldlines is isotropic, H is precisely the Hubble rate. In the homogeneous frw 
case, H is constant at each instant of time t — to- But in the general inhomogeneous 
case, H varies with radial distance from r = on {t = t }. By (|3]), the central 
behaviour of H is 

1 B(w,0)~ 



H(w,y) 



+ 0(y) 



_A(w,0)B(w } 0)_ 

At any given instant w = u>o along {y = 0}, the first term is just the Hubble constant 
H Q = H(wq, 0) measured by the central observer: 

The covariant four-velocity of the dust is given by the gradient of proper time along 
the matter worldlines: = —t,^. It is also given by (H) and (||) as 

«,u = 9^uU v = -Aw,^ +By, fl 
Comparing these two forms, we get 

dt dt 

dt = Adw -Bdy A = —— , B = -— (13) 

ow ay 

which shows that the surfaces of simultaneity for the observer are given in observational 
coordinates by Adw = Bdy. The integrability condition of ([13|) is 



A' + B = (14) 

where a prime denotes d/dy. This is precisely the momentum conservation equation 
(see §). 

Using (D), (|), (0), (HI) and ([14]), it follows that, exactly as in the frw case, the 
Hubble constant is determined by the initial slope of the (e^, z) relation, or equivalently 
of the (dL,z) relation: 

^ -(£)„-(*). < 15 > 



Since the acceleration and vorticity of the dust are zero, the only non-trivial kine- 
matic quantities are the expansion (O) and the shear a. Using (|2|), ( O ) and (|i4|) in 
the general expressions given in H, we find 



a 



V3A 



B 
B 



C 
C 



(16) 



Thus the frw case is given by B/B = C/C. 



3. Reduction of the field equations 

The conservation equations are ([14]) and the energy conservation equation ii M V (U (p) + 
3Hp = 0, which gives, by (0) and (0) 



P 



mN' 
~BCP 



(17) 



The central behaviour of the density is p(w, y) = p(w, 0) + 0(y), and then (0) and (|3]) 
imply 

' ; '+0(y) 



p(w,y) 
N(y) = 



PoBf, 



PoBl 



3m 



B(w,0) s 

y 3 + 0(y 4 ) 



In observational coordinates, the field equations lead to the following first integrals 
(see 0): 



c c 

A + B 

CC A C' 2 AB 
C~C + 2B~C 2 ~ ~ 2C 2 



F = 



N' 



N' 

mN* 
1>CP 



(19) 
(20) 



where N*(y) is an arbitrary function, whose central behaviour is the same as that of 
the number counts: 

~PoB$~ 



N*(y) 



3m 



y" + 0( y 4 



21 



by ©, © and (0). 

In [11, an integration scheme for the field equations, with observational data on 
the past light cone {w = wq}, is presented, intending to give a formal exact solution. 
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Due to the over-restrictive nature of the condition B = A, it fails to do so. Here we 
give a somewhat different approach, which parallels the explicit exact solution in 3 + 1 
coordinates, and which corrects ||. In terms of the intrinsically defined variables [see 



u = it"V M (C) = - , v = (u v k v )- l kfV ll (C) = — (22) 
(II) and © become 

u + v = F 

2 o - m N* 

v + zuv 



C 

Solving and eliminating u, v, we obtain 



A = m (23) 

B = — m (24) 

F-[F 2 -l + mNjC] 1/2 



Thus the problem is reduced to determining C. Integrating (|23|) along the matter flow 
y = const, and using (0), we get 



r 

t + T{y)= m (25) 

where T is arbitrary. This is the ltb exact solution ||, provided we identify 

F 2 = l-kf, k = 0,±l (26) 

Of course, we could have started from this known solution, but then we would have to 
evaluate some of the field equations in observational coordinates anyway in order to 
determine A and B. Furthermore, the observational coordinates allow us to incorporate 
observations directly into the field equations. 

By (E§), is a modified number count variable: 



o 



NM = I y/l - kf(y)*N'(y)dy (27) 

so that k = => N* = N. 
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Now the cosmic time t in (25) is not known explicitly in observational coordinates, 
so we must transform (|25|) into an equation involving only observational variables. This 
entails losing the explicit exact 3 + 1 general solution, but the balancing advantage, 
absent in the 3 + 1 case, is that the light rays are known explicitly, and the fundamental 
equation involves only observational variables. 

We evaluate the integral in (p5|) , differentiate down the light rays, and use ([13]), 



and (|27|) , to get: 



+ 



I mN 3 
C 



C 3 

mN 



T = , for k = 



(2f 



a 



VT^kp 



C 



r 



[mJV,(r-sinr)// 



31' 



T = 2aicsm(fy/C/mN* 
[mN*(smhT - T) / f 3 ]' 



T = 2arcsinh(/ JC/mN, 



where 
for k = 1 

where 
for k = — 



(29) 



For each k, this is the fundamental equation, representing the observational reduction 
of the field equations, analogous to the 3 + 1 case. It is a first-order ode for the area 
distance C as a function of y on each light cone. Though the arbitrary function T'(y) 
generally leads to an infinite family of solutions, on our past light cone {w = wq} 
this fundamental equation relates T'(y) to observationally determined quantities. The 
observational data on {w = Wq} may be taken as cIa(z) and N(z). In summary, we can 
determine the arbitrary functions f(y), N*(y) and T'(y) from the data in the following 
way, very similar to what is suggested in ||: 

(a) Set B = A on our past light cone. This amounts to a choice of y on {w = Wq}, 
and is valid since we do not specify B = A for all w. (If we do not use this freedom, 
then it seems very difficult, if not impossible, to use all the available data on our past 
light cone in determining the solution.) 

(b) Then integrate the Ehlers (or null Raychaudhuri) equation on our past light cone, 
as outlined in || [equations (31) to (33)]. This enables us to change our data from 
functions of z (the observed form) to functions of y on our past light cone. 

(c) We can also determine C on our past light cone, as outlined in |§ (Section 5 iii). 
C is essentially given on our past light cone by the data. Thus, since we know both A 
and B on our past light cone, we can determine f(y) from fli~9|) . 
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(d) Knowing f(y), we can use (f27|) to determine N*(y) since we know N(z) from data 
and thus N(y) from (b). Thus we have determined the arbitrary functions f(y) and 
N*{y) by using the data we presume we have available on our past light cone. 

(e) Finally, as indicated above, we can use (p8|) , (p9|) on {w = w } to find T'(y). 

4. Characterisation of homogeneity 

Observational characterisations of homogeneity (as opposed to unverifiable philosophi- 
cal assumptions such as the Copernican principle ||), are complicated by the fact that 
even the homogeneous frw universes appear inhomogeneous - because observations are 
made down the light cone, which is a hypersurface of inhomogeneity |§ , f| , f| • Homo- 
geneity of an frw universe is implied indirectly by the specific forms of the observed 
quantities. Any deviation from these forms indicates real inhomogeneity (i.e. on spa- 
tial hypersurfaces) . A precise observational characterisation of homogeneity within the 
ltb class of universes has been given 0,0: 

isotropic dust universes are homogeneous if and only if the (area distance, redshift) 
and (number count, redshift) relations measured by the central observer take exactly 
the frw form. 

This characterisation may be slightly weakened in the parabolic (k = 0) case. In 
the fundamental equation fl28|), we may transform the independent variable to N: 

1 - ^mN/C 

It is well known H] that T = const, reduces the metric to FRW. The first order ODE 
(|30|) then shows that d^{N) = C(wq,N) takes the frw form (see section 5 for the 
explicit form) if and only if T = const. Thus, for k = 0, one does not need the (area 
distance, redshift) and (number count, redshift) relations separately to take the frw 
form, but only the single (area distance, number count) relation: 

a parabolic LTB universe is Einstein-de Sitter if and only if the (area distance, number 
count) relation measured by the central observer takes exactly the Einstein-de Sitter 
form. 

This result is consistent with 0, but it follows more directly via our approach. Note 
that the result does not extend to k ^ 0, essentially since in that case 2 arbitrary 
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functions (T and /) need to be fixed by 2 independent observational relations. 

We emphasize the 'global' nature of these characterisations, i.e. the fact that (Ia{z) 
and N(z) must be frw for all z. Local observational parameters derived from d,A 
and ./V are insufficient to characterise homogeneity. For example, the Hubble constant, 
derived from cLa as in (|T5|), cannot in itself distinguish between inhomogeneous and 
homogeneous universes. The dependence of H on the spatial geometry (i.e. on k) and 
on the density p a t the observer, follows from (JT2]), (|2"3"|), and (^) as 



Hn 



\ 3 



P0 k 



f(0) 
fin 



2 



(31) 



The frw case has /'(0) = 1 and B = a (see (|32| ) and (^) below), and for a given p 
and /c, the frw H is observationally indistinguishable from the general ltb H q . 

A particular consequence of ([J]]) arises when an ltb under-dense region is matched 
to an frw region. The low value for po implies a low value for Hq if k = 0, 1. Thus 
cosmological models with a locally high Hubble constant arising from an under-dense 
inhomogeneous local region within a global frw universe (see for example [JTUJ ) , must 
have k = — 1 in the local LTB region. The conditions for smooth matching at the 
boundary (see ]TT[) then require k = — 1 in the frw region, i.e. an open frw universe. 

We now return to the problem of characterising homogeneity, and derive an alterna- 



tive condition (for all k). As pointed out in |]T2|, the assertion in || that the coordinate 
freedom (f|) could be used to set 

B(w,y)=A(w iy ) (32) 

is erroneous, and in fact ( |3"2"D implies restrictions on the spacetime geometry. This choice 
made in || automatically over-restricts the solutions obtained by that procedure to 
be homogeneous (frw) solutions, as we show below. As long as y is specified to be 
comoving, one can set B = A only on one past light cone, i.e. B(wo,y) = A(wo,y), 
and not generally. Even though the ww and wy components of the Lie derivative of 
the metric g^ v along w M are equal, the time derivatives of A and B are not given by the 
Lie derivatives. And there is no further freedom to specify them arbitrarily - as that 
freedom has already been used up by specifying y to be comoving. If, however, y is 
not chosen comoving, then one could set B(w,y) = A(w,y) without loss of generality. 

This correction provides another characterisation of frw universes within the ltb 
class: 
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an ltb universe is frw if and only if there exists a choice of radial comoving coordinate 
such that g wy = g ww in observational coordinates. 

The proof proceeds simply and directly from the formalism set up in sections 2 and 
3. The assumption (|32|) allows us to integrate ([14]) and (|T9|) , to get 



A(w, y) = 0(77) , C(w,y) = a(rj) f \Jl - kf(y) 2 dy , rj = w- 



(33) 



where we have used the central condition (|3|) for C. Then it follows immediately from 
(|33|) that the shear (|T^) vanishes, so that the spacetime is frw. 

For completeness, we present the frw metric and observations in observational 
coordinates. By ( |33"D we can write the metric (0) as 



ds 2 = a{rjY 



-drj 2 + dy z 



l-kf{yfdy\ dn 



showing that 77 is conformal time. It follows that y is conformal proper radial distance 
X, so that [|13|] 



f o V 1 - k f(y) 2d v = sinn 



Now (p2|) , C|33f) and (|34j) reduce the field equation (see [| 



(34) 



C AC 



A 2 ACC" A 2 C" 2 

+ T^^^7 + 



C AC 2C 2 BC C 2B 2 C 2 



to 



which has solutions 



2a- 



d 2 a 
drj 2 



- ha 2 = 



(w - x) 2 k = 

a(w — x) — a { 1 — cos(w — x) k = 1 
cosh(-u; — x) — 1 k — — 1 



(35) 



where a is constant and the big bang occurs at w — \- 

The derivation of the exact observational formulas in frw universes is immediate 
and transparent in observational coordinates. If to is the age of the universe (so that 
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to — Wo, since x = f° r the observer), then the constant a follows from 

as 

1 f 2t ° 3 
a = — < sin t (cos t — 1) 2 

7/0 [ sinht (cosht - 1)~ 2 
while the density parameter O — Po/3Hq follows from ([17]) and (|T8|): 



and fl3"5]) 



(36) 



tt = { 2/(cost + l) 
2/(cosht + 1) 



(37) 



Then (|3~5|) , fl3~6f ) and (|37|) lead directly to all relations between the redshift [see (0)], 
observer area distance [see (|), fl33j)l, luminosity distance [see (§)] and number counts 
[see fllPPI- For example, the number count/ luminosity distance relation is 



32 
3 



1+ l-Vl + 2% 



' 11 fio(l-^o)- 3/2 [sinh2 X K)-2xK)] 



where, for k — ±1: 

' arctan (2v^^T(l + Fo^l) [%(2 - ^oF* 



- arccos fl H d L J^H^dl + (fi - l)(8H d L + 4) 



-1/2 



In 



ftotfcA - 2vT=TWl + 



+ 



In 



Fo^(2 - fio) - 2(1 + H d L )^T=Tk 



(38) 



(39) 



5. Observational generalisation of the Einstein— de 
Sitter universe 

The ltb solution (^5|) allows one to write down any number of explicit exact inhomoge- 
neous generalisations of frw solutions. One simply chooses T(y). However, if we want 
an observationally based generalisation, then the ltb solution is not appropriate, since 
t is not an observational variable, and deriving the observational variables requires in- 
tegration of the null geodesies. In the observational coordinates, the key equation (|28| ) 
or ( p9|) is far more complicated, but it has the advantage that it is already in observa- 
tional variables. Thus any explicit solution will be a readily interpreted observational 
generalisation of an frw universe. 
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We have found one such solution in the case k = 0, i.e. an observational inho- 
mogeneous generalisation of an Einstein-de Sitter universe. The equation ( |30"D can be 
simplified by defining the observational variable 



D 



l_C_ 

mN 



(40) 



Then (|30|) becomes 



3D 

6N 0N +2D+1 



By ©, (HD, (H and (0), we get 



3 dT 
m dN 



l-D 



A(w, N) = 2mND 2 D 

B(w,N) = mMD 2 (D + 2ND N )/(D - 1) 

C(w, N) = mND 2 

p(w, N) = (D- l)/[m 2 N 2 D 6 (D + 2ND N )} 



(41) 



(42) 



Now it is clear that the ODE fl4"T| ) is separable for OT/dN = \mf3, (3 = const. If 
(3 = 0, then T = const., which is the frw case. In the observational variables, the 
Einstein-de Sitter solution following from (f4~TD and (|42| ) is 



A = \mP(P - N 1/3 ) 2 B = ±mMN- 2 / 3 (P - N 1 
C = \mN l l'\P - jVV3)2 p = 192m -2(p _ iv 1 /3)-6 



/3\2 



(43) 



where P(w) is an arbitrary function, subject to P(wq) = 2(mH )~ 1 ^ 3 , which ensures 
that p — 3Pq. The (area distance, number count) relation is cIa{N) = C(wq,N). 
Note that by (0), ([43]) and (|T7|) , the observational variable Z) on {w = wo} determines 
the variation in number counts in terms of the variation in density: 



N p 



The inhomogeneous case T = \m(3N, /3^0, leads to the solution 



(44) 



iV = Q(w) exp 



y/C/mN 



QD 3 dD 



2D 4 + D 3 — (3D + (3 



(45) 
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where Q is an arbitrary function. For small number counts [equivalently, by (|3[), Q2TD 
and (^7|), for large D], the solution fl45 ) on {w = wo} has the form 



d A (N) = {mH^ 2 ) 1/z N 1/z -(m 2 H l )N 2/ * + \mN + 

+ !/3(m 4 tf ) 1/3 iV 4/3 - l(3(m 5 H 2 )N 5 / 3 + 0(N 2 ) 



(46) 



where Q(wq) has been chosen so that the leading term in ([46]) agrees with that in 

(f43|). The first three terms in (^) are exactly the Einstein-de Sitter form given by 

(f43|). Thus by the characterisation proved in section 4, the new solution approaches 

the Einstein-de Sitter solution for small number counts and then deviates from it for 
larger number counts. Choosing y = z on {w = w }, it follows from ([H]) and (p~8|) that 

so that small number counts are equivalent to small redshift. This solution represents 
a universe that is locally homogeneous, but for larger redshift develops real inhomo- 
geneity. If the Copernican principle is not invoked, and if large-scale observations of 
the galactic distribution are taken at face value, then such a model of the late universe 
is not automatically ruled out. 

For completeness, we give the explicit form of the integral in (^), which depends 
on the roots ai of the quartic in the integrand jlJJ: 



(a) (3 = (3± = 13 ± ^rVS «i = a 2 = r = r, a 3 = a 4 = A ^ A ^ r : 



N = Q{w) 1 - 



where 



C 



r V mN 



A 



mN j 



exp 



V 



P 



3r 2 [2r(A + A) - 3AA - r 2 ) 
(A-r) 2 (A-r) 2 

3r 3 [-r(A + A) + AA + r 2 ] 
(A-r) 2 (A-r) 2 

-3A 3 

^~ (A-r) 2 (A-A) 2 



/ C 
mN 



(47) 
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(b) 7^ P 7^ (3± =>■ no ai equal : 



N = Q{w) 1 




(4£ 



where 



and i is not equal to j, k or /. 

Further calculations show that the qualitative behaviour of the solutions (|4"T| ) and 
( pE8| ) falls into two types, according to whether or not the past light cone reconverges. 

6. Concluding remarks 

Motivated by the direct observational evidence for isotropy, and by the need to identify 
indirect observational evidence for homogeneity, we have completed an exact observa- 
tional analysis of isotropic dust cosmologies, i.e. ltb spacetimes. Such an analysis, 
based on the past light cone (which is the hypersurface on which observational data 
is given), was begun in H and continued in ||. An important error in || has been 
corrected here, allowing us to properly complete the integration scheme, as described 
in Section 3. The main result is the reduction of the field equations to a single ODE on 
the past light cone [equations (28|), (p9|)], which then allows us to show in detail how 



the observations, i.e. (area distance, redshift) and (number count, redshift) relations, 
provide the data needed on the past light cone to integrate the field equations and 
determine the metric and matter distribution. 

The observational ODE (p8|) leads to a simplified characterisation of Einstein-de 

Sitter universes via the single (area distance, number count) relation (Section 4). Fur- 
thermore, we give in Section 4 a new observational characterisation of homogeneity, 
i.e. the condition that B = A in observational coordinates. Finally, in Section 5 we 
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used ( P%D to construct an exact observational generalisation of the Einstein-de Sitter 



universe, given by (jig) , (f46|) , (|47|) and (§||). This spacetime is locally homogeneous but 
inhomogeneous on large scales. We are not proposing it as a realistic model of the uni- 
verse, which is more generally believed to be locally inhomogeneous but homogeneous 
on large enough scales, but it is not ruled out by current galactic observations, and the 
method of deriving it could be extended to construct models with different behaviour 
of the inhomogeneity. 
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